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Patterned irradiation of cuprate superconductors with columnar defects allows a new generation 
of experiments which can probe the properties of vortex liquids by confining them to controlled 
geometries. Here we show that an analysis of such experiments that combines an inhomogeneous 
Bose glass scaling theory with the hydrodynamic description of viscous flow of vortex liquids can be 
used to infer the critical behavior near the Bose glass transition. The shear viscosity is predicted to 
diverge as \T — Tbg\~ z at the Bose glass transition, with z ~ 6 the dynamical critical exponent. 



In the mixed state of cuprate superconductors the mag- 
netic field is concentrated in an array of flexible flux 
bundles that, much like ordinary, matter, can form crys- 
talline, liquid and glassy phasesd. The dynamics of the 
flux-line array determines the resistive properties of the 
material and has therefore been the focus of much at- 
tention. Novel types of glasses ase also possible because 
of pinning in disordered samplesn. In particular, the in- 
troduction of columnar damage tracks by heavy-ion ir- 
radiation yields a low-temperature "Bose glass" phase- 
in which every vortex is trapped on a columnar defectEl 
and the rP»ning efficiency of vortex lines is strongly 
cnhancedcTQ. At high temperatures the vortices delocal- 
ize in an entangled flux-line liquid. The high tempera- 
ture liquid transforms into a Bose glass via a second order 
phase transition at Teg, characterized by universal criti- 
cal exponentsEffl. We show here that there are very strong 
divergences in the vortex shear viscosity and other trans- 
port coefficients as this transition is approached from the 
liquid, similar to behavior conjectured for glass transi- 
tions in ordinary forms of matter, and propose experi- 
ments which test our predictions. Vortex matter with 
columnar defects thus provides a concrete example of a 
glassy phase accessed via a genuine second order phase 
transition and characterized by universal critical expo- 
nents. 

The Bose glass transition has been studied theoret- 
ically by viewing the vortex line trajectories as the 
world liflps of two-dimensional quantum mechanical 
particlesad. The thickness of the superconducting sam- 
ple corresponds to the inverse temperature of the ficti- 
cious quantum particles. In thick samples the physics of 
vortex lines pinned by columnar defects becomes equiva- 
lent to the low temperature properties of two-dimensional 
bosons with point disorder. The low temperature phase 
is a Bose glass where the vortices behave like localized 
bosons. It has vanishing linear resistivity and an infinite 
tilt modulus!! The entangled flux liquid phase is resistive 
and corresponds to a boson superfluidcl. 

Although an exact theory of the continuous transition 
at Tbg{B) from the Bose glass to the entangled flux liq- 
uid (or "superfluid" ) is not available, most physical prop- 



erties can be described via a scaling theor y i p . terms of 
just two undetermined critical exponents0'II2lfcy. In the 
low temperature Bose glass each flux line is localized in 
the vicinity of one or more columnar pins. Its excur- 
sion in the direction perpendicular to the applied field 
is characterized by a correlation length that diverges at 
Tbg, l±(T) ~ |T - Tbg\~ v± ■ There is also a diverg- 
ing correlation length along the applied field direction 
(here the^j direction), l^T) ~ \T - T BG \- V W , where 
Uu = 2i>±tB. The time scale r for relaxation of a fluc- 
tuation of size l± is assumed to diverge with a critical 
exponent z, r ~ l z ± ~ \T - T BG \~ ZV± U. The universal 
critical exponents as determined by the most recent sim- 
ulations are v± — 1 and z ~ 4.6 ± 2tL3. Scaling can then 
be used to relate physical quantities to these diverging 
length and time scales. In particular, the resistivity p(T) 
for currents applied in the ab plane is predicted to vanish 
as T -» T BG from above as p ~ \T -T BG \ v ^ z ~ 2) n. Some 
predictions of the scaling theory have been tested exper- 
imentally, but there are as yet no direct measurements 
of the transport coefficients usually associated with glass 
transitions in conventional forms of matter, such as the 
shear viscosity. As we shall see, the behavior of the shear 
viscosity is determined by the dynamical critical expo- 
nent z that controls the divergence of the relaxation time 
in the Bose glass phase. A measurement of the shear 
viscosity would provide a direct probe of the diverging 
relaxation time associated with glassy behavioilij. 

Patterned irradiation of cuprate superconductors with 
columnar defects allows for a new generation of experi- 
ments that may in fact provide a direct probe of^sascous 
critical behavior near the Bose glass transition^. By 
starting with a clean sample, at temperatures such that 
point disorder is negligible, it is possible to selectively 
irradiate regions of controlled geometry. An example 
is shown in Fig. 1. The side regions have been heav- 
ily irradiated, and are characterized by a high matching 
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field and transition curve T B q, while the channel is 

lightly irradiated with a lower matching field B^' < B^ 

and transition curve T B %. When T B % < T BG < T B % the 
flux array in the channel is in the liquid state, while the 
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contacts are in the Bose glass phase. Flow in the re- 
sistive flux liquid region is impeded by the "Bose-glass 
contacts" at the boundaries, as the many trapped vor- 
tices in these regions provide an essentially impenetrable 
barrier for the flowing vortices. As discussed in Ref. [if], 
the pinning at the boundaries propagates into the liq- 
uid channels by a viscous length 5 that depends on the 
flux liquid viscosity. As the temperature is lowered at 
constant field, so that the Bose glass transition T^q of 
the liquid region is approached from above (Fig. 2) the 
growing Bose glass correlations increase 8 and strongly 
suppress the flow in the channel and the associated flux 
flow voltage drop across the channel. 



the local electric field from flux motion at position x in 
a channel of thickness L takes the form 

(1) 

where b > 1 is the length scaling parameter and t = 
\T — Tbg\/Tbg the reduced temperature. This ansatz 
follows from the usual assumption that the continuous 
transition is described by a single diverging length scale 
and the homogeneity condition on the relevant physical 
quantities at the transition (see, e.g., Refs. |^Jll]). 
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FIG. I. A weakly irradiated channel (white region) 
where the flux liquid is sandwiched between two heavily 
irradiated Bose-glass contacts (shaded regions). A cur- 
rent J applied across the channel yields flux motion along 
the channel. The reduced field profiles given by Eq. (9) 
are shown for a few values of S/L and can be measured 
by a series of voltage taps. 

In this paper we analyze experiments with flux flow in 
such confined geometries by combining the predictions 
of the Bose glass scaling theory - generalized to the spa- 
tially inhomogeneous case -.srith the hydrodynamics of 
viscous flow of vortex liquidstJ. Our analysis shows that 
the viscous length 5 controlling boundary pinning is just 
the Bose-glass localization length, and therefore pro- 
vides a prescription for measuring the Bose glass scaling 
near the transition. Both flow in the channel geometry 
sketched in Fig. I and in the Cothino disk geometry (Fig. 
3) used recently by Lopez et alJlj is discussed. Such ex- 
periments can be used to extract the critical behavior 
of various transport coefficients and map out the entire 
critical region. In particular, the flux liquid shear viscos- 
ity is predicted to diverge as \T — Tbg\~ z at the Bose 
glass transition. Because z ~ 4.6 ± 2.0, this powerful 
divergence is reminiscent of the Vogel-Fulcher behavior 
■q ~ exp [c/(T — Tg)] conjectured for glass transitions in 
conventional forms of matter. 

The Bose glass scaling theory summarized earlier is 
easily generalized to the case of spatially inhomogeneous 
flow in constrained geometries. Considering for simplic- 
ity the channel geometry, a generalized scaling ansatz for 
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FIG. 2. A sketch of the (B,T) phase diagram for 
the flux array in the weakly irradiated channel region. 
The heavy line B^ G (T) denotes the continuous transi- 
tion from the Bose glass to the entangled liquid. Also 
shown is the location B B q(T) of the Bose glass tran- 
sition line in the heavily irradiated contacts. When a 
field BgQ (T) < B < BgQ (T) is applied, the flux array 
in the channel is in the liquid state, while the contacts 
are in the Bose glass phase. By decreasing the tempera- 
ture at constant field as indicated by the arrow, the Bose 
glass transition of the channel region is approached from 
above. 

The response in the Bose glass is generally nonlinear 
in the applied current J. By choosing b = t~ v± ~ 1±{T) 
we obtain 

E(T,J,x,L) = l ± (2) 

In the entangled flux liquid the response is linear at small 
current. Upon expanding the right hand side of Eq. (g) 
we obtain for J — > 

E{J^0,x,L) ~ P J— ) 2 ~*JHx/lx>L/lx), (3) 
\a J 

where ao is the vortex spacing and po = (^o^o /cj (l/7o) 
is the Bardeen-Stephen resistivity of noninteracting flux 
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lines, with 70 a bare friction. A scaling form for the re- 
sistivity p{T,L) = AV/(LJ), with AV the net voltage 
drop across the channel, is easily obtained by integrating 
Eq. (H), with the result, 



p(T,L)=p f (T)f(L/l x ) 



(4) 



with f(x) — ~ L du!F(u } x) a scaling function and Pf{T) 
the bulk resistivity, 



V c / 7 



(5) 



In the second line of Eq. (|J) the dependence on the Bose 
glass correlation length l± has been incorporated in a 

renormalized friction coefficient 7 = 70 



(ao) 



that di- 



verges at the transition as 7 ~ \T — Tbg\ vj ~( z ~ 2 ^- For 
L>Ii, the channel geometry has no effect and one must 
recover the bulk result, leading to f(x^> 1) ~ 1. 

The scaling function T can be determined by assum- 
ing that the long wavelength electric field of Eq. (||) 
is described by hydrodynamic equations^. For simple 
geometries where the current is applied in the ab plane 
and the flow is spatially homogeneous in the z direction, 
these reduce to a single equation for the coarse-grained 
flux liquid flow velocity v(r)J13 



- 7v + T^Viv + f L = 0, 



(6) 



The second term in Eq. (^|) is the flux liquid vis- 
cosity rj(T, H) and represents the viscous drag arising 
from intervortex interactions and entanglement. Finally, 
= — ino^oz x J is the Lorentz force density driving 
the flux motion. Intervortex interaction at the Bose-glass 
boundaries translates into a no-slip boundary condition 
for the flux liquid flow velocity. By preventing the free 
flow of flux liquid, the Bose glass boundaries can signif- 
icantly decrease the macroscopic flux-flow resistivity of 
the superconductor. Once the velocity field is obtained 
by solving Eq. (|^) with suitable boundary conditions, 
the electric field profile in the superconductor is found 
immediately from E(r) = ^^-z, x v(r). It is instructive 
to rewrite Eq. (||) as an equation for the local electric 
field, 



<5 2 ViE + E = p f J, 



(7) 



where <5 = yj 77/7 is the viscous length. When the first 
term on the right hand side is absent, i.e., the flux liquid 
viscosity is small, this equation of "viscous electricity" re- 
duces to Ohm's law with flux flow resistivity given by the 
bulk value, pf(T). Interactions, however, make the vis- 
cous drag important and as a result the electrodynamics 
of flux-line liquids is highly nonlocal near the Bose glass 
transition. 

The solution of the hydrodynamic equation for the sim- 
ple channel geometry sketched in Fig. 1, with a homo- 
geneous current J = — xJ applied across the channel, is 
given by 



E(x, L) — pfJ 



1 - 



cosh(x/<5) 
cosh(L/26) 



(8) 



and is shown in Fig. 1. Upon comparing Eq. (g) to Eq. 
(3), we see that the quantity in square brackets in Eq. 
(3) is the scaling function T and find that the viscous 
length 6 is in fact the Bose glass length l±. As the fric- 
tion diverges at Tbg according to 7 ~ \T — Tbg\~ Vj -^ z ~ 2 ' , 
this identification immediately gives that the flux liquid 
shear viscosity also diverges at the Bose glass transition 
with 
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The scaling form for the resistivity is obtained by inte- 
grating Eq. (||), with the result 



p(T,L)=p f (T) 



1 



— — tanh ( — 
L \2l± 



(10) 



If Zj_ <C i, we recover the bulk result of Eq. (|), 
p(T,L) = p f (T) ~ \T- T B g\ u±{z ~ 2) - Near the transi- 
tion, where l±_ L, the resistivity depends on the chan- 
nel width and is controlled by the shear viscosity, with 



P{T,L) 



PfL 2 
121 2 



n (j)o \ 2 L 2 



12rj(T) 



L 2 \T - 



T, 



BG 



(11) 



This strong divergence of the viscosity is precisely the 
kind of behavior expected at a liquid-glass transition. In 
this sense the Bose glass transition is an example of a 
glass transition that is well understood theoretically and 
where precise predictions are available. 
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FIG. 3. Top view of the Corbino disk geometry with 
Bose glass contacts. The magnetic field is out of the page. 
The vortex array is in the Bose glass state in the inner 
and outer densely dotted regions and in the flux liquid 
state in the weakly irradiated annular region. A radial 
driving current drives flux motion in the azimuthal direc- 
tion, and the voltage taps allow the electric field profile 
to be determined. 
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Another important patterned geometry is the Corbino 
disk, recently used by Lopez et al. for defect-free 
materialslia to reduce boundary effects in the flux flow 
measurements. Here we propose fabrication of a Corbino 
disk with Bose glass inner and outer contacts sketched in 
Fig. 3. A current / injected at the outer boundary and 
extracted at the inner boundary yields a radial current 
density J (7') = — 27r (Ri-R 1 ) r that drives vortex motion 
in the azimuthal direction. The electric field induced by 
flux motion is radial, E(r) = —E(r)r, and its magnitude 
is obtained by solving Eq. (Q) in a cylindrical geometry, 
with the result, 



E(r) 



where 



Pfl 



2tt(R 2 - Rx)h 



-+ci/i(f- ) + c 2 Xi(f ) 



K 1 {p 2 )/p 1 - K 1 (p 1 )/p 2 
Kxipx^M - KMhipi) 

h{pi) / P2 - h{p2) / Pi 
^i(pi)/i(p 2 )-^i(p 2 )/i(pi) : 



(12) 



(13) 



Experiments with patterned geometries near the Bose 
glass transformation provide an exciting opportunity to 
probe viscous behavior near a second order glass tran- 
sition. A similar scaling analysis leads to predictions 
for the additional viscosities, which characterize the dy- 
namics of vortex matterli3ii3. For example, the viscous 
generalization of Ohm's law for transport parallel to the 
applied field reads 



(15) 



with pit ~ 1, 



lu , and 5 1 



This work was supported by the National Science 
Foundation at Syracuse through Grants No. DMR97- 
30678 and DMR98-05818 and at Harvard through Grant 
No. DMR97-14725, and by the Harvard Materials Re- 
search Science and Engineering Center through Grant 
No. DMR98-09363. 



with pi.2 = Ri,2/l± and h(x) and K\(x) Bessel func- 
tions. The electric field profiles are shown in Fig. 4. 
The resistivity is defined in terms of the net voltage drop 
AV12 between the inner and outer radii as p{T, R\, R 2 ) = 
AVi 2 /{I/2it(R 2 — i?i)]. Near the Bose glass transition, 
where l± 3> R2, Ri, we find 



P 



4r?(T) 



R\ 



ARlR 2 2 [\n{R 2 /R 1 
K 2 - U 1 



-}. (14) 



As in the channel geometry, the resistivity at the transi- 
tion is completely determined by the diverging viscosity 
and the geometrical parameters of the channel. 




r (fim) 

FIG. 4. The electric field profile for the Corbino disk 
geometry (Eq. (12)), with Ri = 2p,m and R 2 — 12 pm, 
with a width R 2 — Ri = 10pm. The vertical axis 
represents the reduced field 2irE(r)/(pfI) in fim. The 
dashed line is the electric field for vanishing viscosity, 
Eo(r) — 2 7r(fl2~fli) r- sna ded regions represent the 
Bose glass contacts. 
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